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The mammalian cellular microenvironment is shaped by soluble factors and structural components, 
the extracellular matrix, providing physical support, regulating adhesion and signalling. A global, 
quantitative mass spectrometry strategy, combined with bioinformatics data processing, was 
developed to assess proteome differences in the microenvironment of primary human fibroblasts. 
We studied secreted proteins of fibroblasts from normal and pathologically altered skin and their 
post-translational modifications. The influence of collagen VII, an important structural component, 
which is lost in genetic skin fragility, was used as model. Loss of collagen VII had a global impact on 
the cellular microenvironment and was associated with proteome alterations highly relevant for 
disease pathogenesis including decrease in basement membrane components, increase in dermal 
matrix proteins, TGF-p and metalloproteases, but not higher protease activity. The definition of the 
proteome of fibroblast microenvironment and its plasticity in health and disease identified novel 
disease mechanisms and potential targets of intervention. 
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Introduction 

Mammahan cells secrete soluble and structural proteins into 
the extracellular space. The structural components make up 
the extracellular matrix (ECM) and both, matrix and soluble 
proteins, such as growth factors, proteases and cytokines, 
shape the cellular microenvironment and regulate cell 
functions. Transmembrane receptors like integrins connect 
the extracellular space with intracellular signal transduction 
relays. The ECM regulates not only cell adhesion, but also 
bioavailability of growth factors and composition of mem- 
brane receptors (Brizzi et aZ, 2012). The human skin consists 
of an epithelium, the epidermis and a subjacent mesenchymal 
tissue, the dermis, which both generate complex and 
functionally specialised ECMs: the subepidermal basement 
membrane and dermal networks of fibrillar collagens and 
associated proteins. Cell-ECM interactions in the skin have a 
vital role in tissue homeostasis, wound heahng and ageing, 
and are often deregulated in disease (Watt and Fujiwara, 
2011). Hence, its importance and its easy accessibihty make 
skin ECM highly attractive for studying ECM composition 
and plasticity. 



Physiological skin integrity and resistance to mechanical 
stress rely on the functions of the dermal-epidermal junction 
zone (DEJZ), which anchors the epidermis to the underlying 
dermis. Protein complexes at the DEJZ, such as hemidesmo- 
somes and focal adhesion complexes, mediate interactions of 
the cytoskeleton in basal keratinocytes with the basement 
membrane. Anchoring fibrils, which emanate from the base- 
ment membrane into the dermis and entrap dermal collagen 
bundles, estabhsh stable dermal-epidermal adhesion. The 
major component of the anchoring fibrils is collagen VII (C?) 
(Has and Bruckner-Tuderman, 2006} . 

Loss of C7 causes recessive dystrophic epidermolysis 
bullosa (RDEB), an inherited monogenic skin disorder 
with diminished dermal-epidermal adhesion, mechanically 
induced skin bhstering and subsequent scarring (Fine et aU 
2008). Mutations in the C7 gene, C0L7A1, are associated with 
a broad range of clinical phenotypes, ranging from locahsed 
manifestations to generahsed mucocutaneous bhstering and 
scarring, and to high risk of aggressive skin cancer (Kern et aU 
2006; van den Akker et aU 2011). Genotype-phenotype 
correlation studies on RDEB have shown that null mutations 
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and absence of C7 are associated with severe RDEB, whereas 
hypomorphic mutations cause milder RDEB forms 
(Hammami-Hauash et aU 1998). Intriguingly, there are 
significant interindividual differences in patients with iden- 
tical C0L7A1 mutations (Kern et al 2009). 

Very little information exists on the consequences of loss of 
C7 at the cellular level and in relation to the cellular 
microenvironment. Loss of the structural function of C7 
perturbs its interaction with laminin-332, which provides 
dermal-epidermal anchorage and is also required for kerati- 
nocyte survival (Waterman et aU 2007} . In vivo, C7 deficiency 
was associated with increased TGF-pi and accumulation of 
dermal ECM proteins in skin (Fritsch et aZ, 2008). In three- 
dimensional cultures in vitro, C7-deficient keratinocytes 
became invasive, and recent studies suggested that 
RDEB fibroblast ECM, like cancer-associated fibroblast 
ECM, drives development of cutaneous squamous cell 
carcinoma (Martins et al, 2009; Ng et al, 2012). Additionally, 
differences in collagen breakdown by matrix metallopro- 
teases (MMPs) may contribute to distinct RDEB phenotypes 
(Titeux et al, 2008). 

Proteomic analysis using mass spectrometry (MS) repre- 
sents a state-of-the-art unbiased and global approach 
(Zimmermann et al, 2010; Engelke et al, 2012) to questions 
relating to disease pathogenesis. Few studies have investigated 
secreted proteins of skin cells, including human dermal 
fibroblasts, mouse keratinocytes and human oral squamous 
cell carcinoma lines (Boraldi et al, 2003; Pflieger et al, 2006; 
Del Galdo et al, 2010; Todorovic et al, 2010), but the analyses 
were mainly of qualitative nature and the number of identified 
proteins low. In a recent systems biology study, the role of 
MMP2 in skin inflammation was investigated by employing 
isobaric tags for relative and absolute quantitation (iTRAQ) 
highhghting the power of global quantitative proteomics (auf 
dem Keller et al, 2013). Here, we developed a comprehensive 
approach by employing stable isotope labelling by amino acids 
in cell culture (SILAC) -based quantitative MS (Ong etal, 2002), 
followed by bioinformatics analyses, to investigate quantita- 
tive differences of secreted proteins and respective post- 
translational modifications (PTMs) in human dermal 
fibroblasts derived from normal and pathologically altered 
skin, C7-negative RDEB being our model system. We identified 
and quantified patient-specific and general alterations of the 
extracellular proteome upon C7 loss and generated an atlas of 
oxidative PTMs important for structural stability, providing 
a valuable resource of proteins shaping the cellular micro- 
environment. Hitherto unknown proteins, molecular path- 
ways and mechanisms involved in RDEB pathogenesis were 
identified. 

Results 

Secretome of normal skin fibroblasts 

For proteomic profihng of proteins secreted by primary 
human skin fibroblasts, cells were processed as shown in 
Figure lA. Secreted proteins were divided into two fractions: 
(a) cell conditioned medium (CM) comprising all soluble 
proteins and (b) insoluble proteins remaining on the culture 
plate after cell lysis. The latter was defined as ECM in the 
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experimental setting used here. The CM was harvested at day 
4, and the ECM was purified at day 6 by removing the cells 
with basic Triton-NH40H (Vlodavsky, 2001), followed by 
protein extraction with 4% SDS. Compared with other 
isolation procedures this approach yielded the best results, 
i.e., a fibrillar matrix visible in the hght microscope after 
cellular detachment, and enrichment of extracellular proteins, 
as judged by the respective GO (Supplementary Figure SI). 
The efficiency and purity of the ECM was confirmed by 
western blot (Figure IB). 

To globally characterise the composition of the fibroblast 
secretome, the proteins in the ECM and CM were fractionated 
by SDS-PAGE, in-gel digested with trypsin and analysed by 
LC-MS/MS, resulting in the identification of almost 2000 
proteins with a false discovery rate (FDR) below 1%. 
Since the sensitivity of modern mass spectrometers is 
extremely high, it is likely that small amounts of contaminat- 
ing intracellular proteins would be detected in the biochemi- 
cally enriched fraction of extracellular proteins. To 
distinguish these and to generate a comprehensive Ust of 
ECM proteins, we employed a bioinformatics data processing 
pipeline using the following parameters (Figure IC): (1) 
proteins carrying GO terms 'extracellular' (extracellular 
region, space and matrix) were filtered; (2) to identify new 
ECM components, all proteins reproducibly identified in the 
ECM and the CM were also considered as extracellular; 
(3) proteins predicted by in silico analysis to carry a 
signal peptide (SignalP; www.cbs.dtu.dk/services/SignalP/) 
(Petersen et al, 2011) were counted as extracellular 
(Henningsen et al, 2010). These measures yielded ^660 
extracellular proteins in the ECM and 740 in the CM from the 
total list of identified proteins, with >60% being detected in 
both fractions (Figure IC). The filtered proteins were 
analysed based on their Swiss-Prot (SP) and Protein 
Information Resource (PIR) keywords (Figure ID). While 
membrane proteins were enriched in the ECM fraction, 
proteins with enzymatic activities, such as proteases and 
hydrolases, were enriched in the CM, indicating the different 
nature of the two compartments. As expected, the terms 
secreted, ECM, and signal peptide, among others, were 
common to both groups. 

To assess potential differences in the abundance of 
extracellular proteins in a physiological setting, we SILAC 
labelled skin fibroblasts of three healthy donors (Sprenger 
et al, 2010). CM and ECM were purified and the data 
processed as outlined above. In two biological rephcates of 
ECM and CM, respectively, we quantified 863 potential 
extracellular proteins, of which 40 % were annotated as being 
extracellular based on GO terms. We observed only minor 
donor-specific differences in ECM and CM, indicating that 
three samples were sufficient to capture proteome alterations 
in the used experimental settings (Figure IE; Supplementary 
Tables SI and S2) . In all, 95 % of proteins were in the interval 
of ± 0.75 (log2 SILAC ratios), not showing altered abundance 
in the different samples, and biological replicates 
showed good reproducibility (r = 0.6-0.93; Supplementary 
Figures S2 and S3). Thus, we established an experimental 
platform for global analysis of protein abundance differences 
of extracellular proteins, ECM and CM, and showed 
that extracellular proteins of primary skin fibroblasts of 
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Figure 1 Experimental (design and data processing pipeline. (A) Primary skin fibroblasts were kept in culture for 6 days. At days 2, 3, and 4, fresh culture medium 
containing ascorbate was added. The CM was harvested at day 4 and the ECM at day 6. (B) Intracellular proteins representing the cell compartments cytosol (Tubulin p 
chain, GAPDH), mitochondrion (GRP-75) and golgi apparatus (golgin-97) are only detected in the cell lysate (i). Procollagen I was in the cell lysate, whereas mature 
collagen I was found only in the ECM (e). (C) The ECM and CM were analysed by MS-based proteomics, and secreted proteins were filtered on (I) GO annotation 
'extracellular', (II) proteins that were repeatedly detected in the ECM and the CM, and (III) proteins detected only in the ECM or the CM if they were predicted to have 
signal peptides. (D) SP and PIR keywords were used to highlight enriched protein classes (BH-adjusted P-value <0.05). (E) Experimental variability was analysed by 
quantitatively comparing protein abundances of differentially SILAC-labelled ECM and CM from the same donor (ctrl 1 and Ctrl 3 in grey, respectively). Biological 
variability was analysed by comparing protein abundances of the ECM and the CM from different donors (in green). Figures were produced using Servier Medical Art 
(www.servier.com). 



healthy donors exhibit only minor differences in their 
secretion pattern. 

Interindividual variance of secreted proteins in 
fibroblasts upon loss of C7 

To study the effects of loss of C7 on the composition of the 
secretome, skin fibroblasts of four individuals with RDEB and 
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three healthy donors were SILAC labelled, subcultured to the 
same number of passages and subjected to the quantitative 
proteome analysis as described above. The diagnosis of RDEB 
was estabhshed by two criteria: complete lack of C7 in the skin 
and C0L7A1 mutations leading to a premature stop codon 
(Supplementary Figure S4; Table I). The ECMs of the three 
controls were combined to generate a Super-SILAC mix, 
minimising the interindividual influences of the healthy 
donors (Geiger et al, 2010} . This mix was then spiked in equal 
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Table I Characteristics of RDEB patients in this study 

Patient Age at Sampling C0L7A1 Consequence Reference 
biopsy location mutation 
(days)/ 
gender 



1. 4/f Thigh 425A>G PTC Kern a/ 

682+1 PTC (2009) 
G>A 

2. 36/m Thigh 425A>G PTC Kern eta/ 

c.5261dup p.G1755RfsX17 (2009) 

3. 6/f Undisclosed 1732C>T R578X Kern eta/ 

1732C>T R578X (2006) 

4. 8/m Undisclosed cl934del p.P645QfsX45 Kern eta/ 

cl934del p.P645QfsX45 (2009) 



amounts into the ECMs purified from medium and heavy 
labelled RDEB cells (Figure 2A). Furthermore, Super-SILAC 
samples were also generated of ECMs labelled control cells 
using the same procedure as for the patients. Subsequently, 
samples and data were processed as outlined (Figures 2A, B 
and IC] . The same workflow was performed for the CM. Since 
the Super-SILAC mix was used as a common standard, it was 
possible to directly compare quantitative differences between 
the four RDEB and the three control samples. We quantified 
587 potential extracellular proteins (45% carrying 'extracel- 
lular' GO terms) in a total set of 190 LC-MS/MS analyses 
comprising at least two biological rephcates for all conditions 
(Supplementary Tables S3 and S4; Supplementary Figures S5 
and S6). Of these, 154 proteins were identified only in CM 
samples (Supplementary Table S4). On average, 45.7% of the 
ECM proteins and 31.7% of the CM proteins of RDEB 
fibroblasts showed abundance differences larger than 0.75 
(log2 SILAC ratio), much larger proportions than the 5% 
observed for the control fibroblasts. 

As all four RDEB patients carried null mutations in the 
C0L7A1 gene, the question was raised whether the ECMs of 
RDEB fibroblasts generally show the same trend of alterations 
in their proteomic composition. To address this, we compared 
protein abundance differences by performing a principal 
component analysis (Figure 2C). Control and RDEB samples 
clearly separated, indicating general alterations in ECM 
composition in RDEB. However, it became evident that RDEB 
ECMs were more heterogeneous than control ECMs and 
showed patient-specific effects. 



Modulated basement membrane and ECM 
proteomes in RDEB 

Quantitative differences of secreted proteins in C7-deficient 
RDEB fibroblasts were analysed by hierarchical clustering of 
ECM (Figure 2D) and CM data (Supplementary Figure S7A). 
Generally, proteins carrying GO terms 'extracellular structure' 
and 'matrix organisation' were enriched in RDEB ECMs, 
whereas proteins related to 'cell morphogenesis' and 'cytos- 
keleton organisation' were less abundant. We examined 
proteins carrying the GO terms 'cell adhesion' and 'extra- 
cellular structure and organisation' in detail. For example, the 
interstitial collagens were increased in the ECMs of patients 1, 
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3 and 4. In contrast, in samples of patient 2, the polypeptides 
encoded by COLlAl, C0L1A2, C0L3A1, C0L5A1, COL5A2, 
C0L12A1 and C0L15A1 did not accumulate. Also the collagen 
modifying enzymes lysyl oxidase homologs 1, 2, and lysyl 
hydroxylase 1, as well as the elastin microfibril interface- 
located proteins 1 and 2, were less abundant in the samples of 
patient 2 as compared with the other three patients 
(Supplementary Figure S8). The complete data sets demon- 
strating patient-specific differences are presented in 
Supplementary Tables S3 and S4. 

Clustering of protein abundances in CM was in agreement 
with ECM findings separating RDEB and control samples 
(Supplementary Figure S7A). However, generally, the ECM 
values spread more than the CM values, and the respective 
data sets correlated only weakly (r=0.25). This could reflect 
altered protein-protein interactions and integration of indivi- 
dual proteins into the ECM suprastructures in RDEB samples 
(Supplementary Figure S7B). Highly abundant proteins in 
RDEB CM were TGF-pi, and the TGF-(3-induced protein ig-h3. 

Quantitative proteomics-based findings were validated with 
several methods. Western blot analyses of RDEB ECMs 
revealed less deposition of basement membrane proteins 
laminin pi, laminin yl, nidogen 1 and fibulin 1, whereas 
tenascin C and TGF-pi -induced protein ig-h3 were accumu- 
lated (Figure 3 A) . These differences correlated well with SILAC 
MS ratios (r = 0.83; Figure 3B). Also immunofluorescence 
stainings of ECM were in agreement with proteomics 
data (Figure 3C; Supplementary Figure S9C). Interestingly, 
laminin (31 and nidogen 1 were less abundant in ECMs but not 
in respective RDEB cell lysates, suggesting altered protein 
trafficking/ECM integration (Figure 3A and C; Supplementary 
Figure S9A-C) . The corresponding in vivo data were congruent 
with the in vitro results. Immunofluorescence staining showed 
increased tenascin C and decreased collagen IV and laminin 
Pl/yl deposition at the DEJZ in RDEB skin (Figure 3D). 
Paraffin-embedded specimens of normal and RDEB skin were 
assessed regarding the pattern and the intensity of Elastica van 
Gieson staining of the dermis. In control skin, the dermal 
ECM presented with a loose, undulating fibrillar network, 
incorporating regular interconnected elastic fibres (Figure 3E) . 
In RDEB skin, the ECM was tightly packed, and the 
fibre orientation paralleled the skin surface, indicating 
quantitative increase in ECM proteins, as detected by the 
proteomics approach, and perturbed organisation, as seen in 
fibrosis. Quantitation of the normahsed mean dermal staining 
intensity revealed a significant fibrotic change in RDEB dermis, 
as compared with locahsation and age-matched control 
dermis. 

To test whether the differences in protein abundance were 
directly linked to loss of C7, siRNA was used to knock down C7 
in normal primary skin fibroblasts. This resulted in a decrease 
in laminin yl and an increase in collagen I in respective cell 
lysates as compared with controls and in agreement with 
proteomics findings (Supplementary Figure S9D and E; 
Supplementary Table S3). 

Finally, linear models for microarray data (LIMMA) (Smyth, 
2004) were used to identify consistent changes in ECM 
compositions between all RDEB and control fibroblasts and 
identified 214 proteins as significantly altered in all RDEB 
ECMs (P<0.05, BH adjusted; Supplementary Table S3). We 
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Figure 2 Analysis of quantitative differences in ECM compositions due to loss of C7. (A) Controls (light) were combined to a Super-SILAC mix to minimise differences in 
reference ECM isolations. These were spiked in 1 :1 :1 ratio to SILAC-labelled ECMs of RDEB fibroblasts (medium and heavy). The samples were processed as outlined and the 
resulting peptide mixtures were analysed by LC-MS/MS. (B) Peptide signal intensities recorded in MS spectra allowed relative quantification of RDEB and control ECM peptides, 
which were then combined to protein ratios. Coloured circles represent the respective SILAC state, as indicated in (A) (exemplified data). (C) Log2-transformed abundance 
changes of filtered ECM proteins were used for a principal component analysis. The two experimental groups are clearly separated, RDEB samples showing a bigger diversity in 
component 2. (D) Protein ratios were log2 transformed and z-score normalised. Columns containing data from the different samples were hierarchically clustered and rows 
containing protein entries were clustered by k-means. The proteins in k-means cluster were tested by DAVID (Huang da etal, 2009) for enrichment of GO CC and GO BP terms 
and selected terms are shown (full list in Supplementary Table S3; BH-adjusted P-value <0.01). Figures were produced using Servier Medical Art (www.servier.com). 



performed a network analysis to assess the functional 
interactions of these extracellular proteins. This generated a 
protein interaction network consisting of 64 members 
(Figure 4). It became evident that basement membrane 
proteins were significantly less abundant in RDEB (P<0.02, 
BH adjusted), in particular, the laminin chains oc4, pi, P2 and 
yl, nidogen 1, collagen IV and several integrins. In contrast, 
interstitial collagens III, V and VI, protease inhibitors like 
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serpins El and E2, and metalloproteases were more abundant 
(P< 0.008, BH adjusted). 

Altered abundance of PTMs of secreted proteins 

We scanned secreted proteins for the presence of hydroxypro- 
line and hydroxylysine residues, which are important for ECM 
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Figure 3 Altered deposition of EClVI proteins in RDEB. (A) Purified EClVI of control (ctrl) and RDEB (pt) fibroblasts were analysed by western blot using antibodies to 
proteins, which were indicated by MS-based proteomics to be significantly altered. Samples were normalised to cell number. (B) Correlation of western blot and SILAC 
ECM protein fold changes of RDEB versus control samples. (C) ECM of control (ctrl) and RDEB (pt) fibroblasts were stained using antibodies against annotated proteins, 
which were indicated by MS-based proteomics to be significantly altered. Fibronectin staining was used as a control for the ECM. (D) Immunofluorescence staining of 
control (Ctrl) and RDEB (pt) skin in the skin of patient 3. (E) Fibrotic dermal ECM in RDEB. Normal skin (left panel): dermal ECM (red) with a loose, undulating fibrillar 
network, incorporating regular interconnected elastic fibres (black), vessels and adnex structures. RDEB skin (middle panel): tightly packed dermal ECM, fibre orientation 
parallel to the skin surface. Lower panels: magnification of detail (white frame in upper panels). Quantitation of normalised mean dermal staining intensity (right panel)— 
as a measure for fibrotic change— (normalisation to epidermal staining, y axis: relative units) shows a significant difference between normal and RDEB dermis (P<0.01 , 
two-tailed, unpaired T-test). Box-plot of results from skin specimens of RDEB patients and age- and localisation-matched healthy controls. n = 5. Scale bars represent 
100|im (C) and 50|im (D). 



stability and cross-linking (Monnier et aU 2005; Loenarz and 
Schofield, 2011), and analysed surrounding amino-acid motifs 
to gain further insights into potential mechanisms and 
biological functions of respective source proteins (Figure 5; 
Supplementary Figure SIO). In total, 634 hydroxyproline and 
130 hydroxylysine sites were identified (locahsation score 
^0.75, class I site, Andromeda score ^100} (Olsen et aU 
2006} . As expected, the largest groups of hydoxyproline and - 
lysine sites were variations of the collagen motifs Gly-X-Pro- 
OH and Gly-X-Lys-OH. However, also new hydroxylation 
motifs in other protein classes as defined by the matrisome 
project (Naba et aU 2012) were identified, especially in the case 
of proteins present in CM, of which many were glycoproteins, 
e.g., fibrillin-1 and -2 (Figure 5 A and B; Supplementary Figure 
SIO). Collagens detected in ECM and CM differed in their PTM 
status depending on the protein and not on the source 
compartment, suggesting that PTMs might have a role in 
protein-matrix integration and are dynamically modulated 

6 Molecular Systems Biology 2013 



(Figure 5C; Supplementary Figure SIOC). Interestingly, hydro- 
xyproline sites identified in ECM and CM proteins in RDEB 
samples differed substantially from respective controls. 
Whereas sites detected in samples from patients 1, 3 and 4 
were on average 2. 3 -fold less abundant than in controls, 
samples from patient 2 exhibited an even higher overall 
abundance level, however with differing contributions of 
single sites (Figure 5D; Supplementary Figure SIOD; 
Supplementary Tables S5 and S6). Thus, loss of C7 in RDEB 
affects not only extracellular proteome composition but also 
the PTM status of secreted proteins. 



Inactivation of metal loproteases 

The protein interaction network shown in Figure 4 indicates 
high abundance of MMP2 and MMP14 and protease inhibitors, 
like TIMP3, in ECMs of RDEB skin fibroblasts. Increased 
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Figure 4 Network analysis reveals decrease in basement membrane proteins and increase in interstitial collagens. Network nodes represent quantified ECM proteins 
and edges known protein-protein interactions. Proteins with significantly lower abundance in RDEB than in controls are coloured blue and proteins with higher 
abundance yellow. Grey proteins show only significant changes in some of the patients. The network was generated by Path Designer, Ingenuity Systems. 



abundance of both proteases in RDEB ECM was validated by 
western blot (Figure 6A) . To estimate the proteolytic activity, 
control and RDEB fibroblasts were seeded on 3D collagen I gels 
and incubated for 3 days. During this period, control 
fibroblasts degraded the gel, whereas RDEB cells showed 
strongly reduced activity. Even after TNFa stimulation, which 
increases MMP14 expression and thereby MMP2 activation 
(Han et aU 2001), only little degradation of collagen I by RDEB 
fibroblasts was observed (Figure 6B). Protease inhibitor 
profiling linked the proteolytic activity to metalloproteases. 
In control cells, TAPI-0 inhibited collagen I degradation, but 
neither leupeptin, nor the serine protease inhibitors SBTI and 
PMSF had an effect (Figure 6B}. To specifically block the 
activity of membrane-bound metalloproteases, the furin 
inhibitor D-RVKR-CMK was included (Overall and Lopez- 
Otin, 2002). D-RVKR-CMK inhibited collagen I degradation by 
control cells, indicating that the proteolysis was due to MMP14 
activity. 

To test whether TIMP3, which can block MMP14 activity 
(Pavloff et al 1992; Will et al, 1996) and accumulated in RDEB 
samples, was involved in MMP14 inactivation recombinant 
TIMP3 was also added to control cells and indeed blocked 
collagen I degradation (Figure 6B}. On the other hand, the 
addition of RDEB CM, which contained soluble protease 
inhibitors, like TIMPl and TIMP2, to control cells had no effect 
(Supplementary Figure Sll). 

The variation of MMP14 activity and abundance may also 
reflect differences in MMP14 recycling (Maquoi et aU 2003} . To 
address the influence of lysosomal degradation of MMP14, 



cells were incubated with concanamycin A (ConA), an 
inhibitor of the lysosomal H^-ATPase. In contrast to RDEB, 
ConA led to a pronounced accumulation of MMP14 in control 
cells (Figure 6C), which was reflected by an increase in 
MMP14 and MMP2 activity (Figure 6B and C), indicating 
impaired MMP14 turnover in RDEB. 

To analyse proteolytic processing of membrane-bound 
MMP14, we fractionated two control and two RDEB ECMs by 
SDS-PAGE. The gel lanes were cut into pieces, ordered 
according to the molecular weight, in-gel digested by trypsin, 
and analysed by MS (Dix et aU 2008). In control samples, the 
MMP14 protein was almost entirely present in the full-length 
form (Figure 6D}. In contrast, in RDEB samples, the full-length 
form represented only about half of the MMP14, and several 
smaller molecular weight forms were observed. Mapping of 
the peptides revealed that the smaller forms consisted only 
of the hemopexin-like domains 1-4, the transmembrane 
region and the cytoplasmic tail, i.e., had undergone shedding 
and lost the catalytic domain (Figure 6E}. Thus, increased 
MMP14 shedding, increased abundance of inhibitors like 
TIMP3, and altered protein turnover contribute to diminished 
MMP14 activity in RDEB cells (Figure 6F}. 



Discussion 

The ECM forms a complex meshwork around cells and organs, 
which does not act as a mere structural support but actively 
shapes the cellular microenvironment and has a vital part in 
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intracellular signal initiation and transduction (Watt and 
Fujiwara, 2011). Here, we employed MS-based proteomics 
relying on SILAC labelling to investigate the protein micro- 
environment and the PTM status of secreted skin fibroblast 
proteins in health and disease. Several hundred proteins and 
their oxidative PTMs, which are important for structural 



stability, were quantified in ECM and CM samples of primary 
dermal fibroblasts. 

Severe RDEB, associated with complete loss of C7, was 
chosen as a model disease to study the influence of a single 
protein on cellular microenvironment plasticity. Since the 
cellular microenvironment has a complex structure and 
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Figure 5 Alterations of oxidative PTIVIs. Hydroxyproline (A) and liydroxylysine (B) sites were scanned for sequence motifs using IVIotif-X (Scliwartz etal, 2009; Ciiou 
and Scliwartz, 2011). Site numbers and respective source proteins grouped according to the matrisome project (Naba et al, 2012) are shown as bar diagrams. 
(C) Numbers of hydroxyproline sites detected in collagens in ECM and CM, respectively, are indicated as bars. (D) Hydoxyproline sites ratios were normalised to 
respective source protein ratios (log2, z-transformed). Sites detected in RDEB and control ECMs are depicted as heatmap. Samples were hierarchically clustered and 
sites by k-means. 



Figure 6 Decrease in protease activity in C7-deficient fibroblasts. (A) Western blot of MMP2 and MMP14 of ECM. (B) The blue rounds represent the intact collagen I 
matrix. Upon degradation of the matrix, the staining becomes white. In contrast to the RDEB cells (pt), control cells (ctrl) degrade the matrix, and the degradation is 
stimulated by the addition of TNFa. Inhibitor profiling demonstrates that the degradation is mediated by metalloproteases, since TAPIO, D-RVKR-CMK and TIMP3, but 
not serine and cystein protease inhibitors leupeptin, SBTI, and PSMF prevent it. ConA, inhibiting lysosomal degradation, leads only in control cells to enhanced 
degradation. (C) ConA leads to enhanced MMP2 activity in control cells as analysed by zymography. Western blot analysis reveals impaired MMP1 4 turnover in patient 
cells compared with control cells. Numbers indicating fold change as compared with untreated control. Samples were normalised to protein amount. (D) ECM proteins of 
RDEB (ptl , pt3) and control (ctrl2, ctrl3) cells were separated by SDS-PAGE and analysed by MS. The total number of detected MMP14 peptides in one sample was set 
to 100% and the distribution of peptides was mapped to respective gel slices. MMP14 appears to be shed in RDEB cells and not in control cells, as only -50% of the 
MMP14 peptides were detected in the 60-kDa gel slice harbouring full-length MMP14 in RDEB cells in contrast to almost 100% in control cells. (E) Schematic 
representation of the MMP14 molecule. Peptides derived from control cells (green) and RDEB cells (yellow) were mapped to the MMP14 molecule. Combined results of 
two biological replicates are shown. (F) Model of altered MMP14 turnover. Under healthy control conditions, C7 is present and active in the extracellular environment. 
MMP14 recycling and degradation is ensured to maintain degradation of collagen I and activation of MMP2. In RDEB, the ECM lacks C7, leading to an accumulation of 
MMP1 4, MMP2 and TIMP3. MMP1 4 and MMP2 activity is impaired due to the presence of inhibitors like TIMP3 and increased shedding of MMP1 4. Additionally, MMP1 4 
recycling is perturbed further amplifying its accumulation. 
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supramolecular organisation, it may not be adequately 
represented in analyses relying on gene expression arrays. 
Therefore, we took an unbiased proteomics approach to 
elucidate molecular disease mechanisms. Compared with a 
recent study using mRNA expression profiling, which detected 
upregulation of C0L5A1, C0L12A1, ITGA3, ITGA6 and TSPl 



genes in RDEB (Ng et al, 2012), we identified substantially 
more deregulated extracellular proteins. Whereas mRNA 
abundances of the respective study correlated weakly with 
intracellular protein abundances of our samples, mRNA levels 
cannot be employed to deduce ECM protein abundances 
(Supplementary Figure SI 2). Hence, the proteomics 
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experiments targeted directly the perturbed cellular micro- 
environment and generated data, which are not accessible by 
transcription analyses. 

So far, only few studies have addressed the complexity of the 
skin fibroblast ECM proteome, often in a mostly qualitative 
manner. We have previously shown that SILAC labelling and 
cultivation of primary skin cells has only little influence on 
their proteomic composition (Sprenger et aU 2010} and we are 
confident that our cell culture system reflects the in vivo 
situation. In the present study, this was corroborated by the 
morphological analyses of normal and RDEB skin. 

The high sensitivity of modern mass spectrometers led to 
detection of intracellular proteins despite careful purification 
of extracellular proteins before proteomic analysis. Therefore, 
a data processing pipeline was designed to filter extracellular 
proteins by bioinformatics means. As our main aim was the 
identification of pathological changes in the ECM proteome, 
we took a non-stringent approach and characterised in total 
872 proteins as being extracellular based on their respective 
GO (39% carrying 'extracellular' GO terms), their repeated 
identification in biological rephcates, and the presence of 
potential signal peptides. Twenty percent of the ECM proteins 
identified in this study were recently defined as components of 
the matrisome (Naba et aZ, 2012). This may not represent a 
final proof that all proteins in our tables are exclusively 
extracellular, but the definition of extracellular proteins is 
indeed not a trivial issue; new extracellular functions have 
been reported for classical intracellular proteins like glucose-6- 
phosphate isomerase, HSP90a and calreticulin (Butler and 
Overall, 2009). 

We identified 764 oxidative PTM sites within the set of 
potential extracellular proteins. PTMs alter protein activity and 
stability, and great effort has been put into their mass 
spectrometric analysis. Enzymatic oxidations of ECM proteins 
are known to be essential for collagen stability, but hydroxyla- 
tions have also been shown to alter signal transduction and 
protein-protein interactions (Loenarz and Schofield, 2011). 
New sequence motifs for hydoxyproline, which differ from the 
canonical collagen motifs and localise to other protein classes, 
could be characterised hinting at so far unstudied molecular 
functions. PTM abundances differed between RDEB and 
control samples, and it is very likely that the alterations exert 
direct effects on ECM stability and structure. 

The disease proteomics approach uncovered new molecular 
players in skin fragility. Many of them have a logical link to the 
major clinical signs of severe RDEB, i.e., weak dermal- 
epidermal adhesion and dermal fibrosis, but the number and 
complexity of contributing molecules is much higher than 
anticipated. It became clear that lack of C7 does not simply 
represent a lost structural link, but causes a number of parallel 
changes on protein and PTM level, which influence the 
dermal-epidermal microenvironment. 

Three novel aspects have clinical imphcations. First, the 
reduced abundance of basement membrane components is 
unexpected and intriguing. Relative deficiency of collagen IV, 
nidogen 1, or the laminin pi and yl chains, which assemble to 
laminin 511, weakens the basement membrane suprastructure 
and amplifies the skin fragility caused by lack of functional 
anchoring fibrils. This is a pivotal finding that will influence 
the design of future molecular therapies. 
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Second, the inherently high production and secretion of 
ECM proteins in RDEB fibroblasts, including tenascin C, and 
interstitial collagens I, III, V and VI, clearly contribute to 
excessive scarring and dermal fibrosis in RDEB. Moreover, 
TGF-(3, a growth factor that stimulates expression of many 
ECM proteins (Akhmetshina et aZ, 2012), was increased in the 
CM of RDEB fibroblasts. The combination of inherently 
elevated ECM synthesis and high secretion of TGF-P by RDEB 
fibroblasts is likely to be a major driver of scarring in RDEB 
skin, presumably further stimulated by other factors, such as 
inflammatory cytokines released through trauma-induced 
separation of the skin layers (Wynn and Ramalingam, 2012). 

Third, the significant increase in proteases and protease 
inhibitors in RDEB ECM re-opens a field of investigation 
relating to molecular disease mechanisms. Skin fibroblasts 
express MMP5 1, 2, 3, 9, 10, 13 and 14 (Loffek etaU 2011) and, 
in the past, investigators have suggested that increased MMP 
activity contributes to RDEB pathology (Winberg et aU 1989). 
The gene MMPl was the first one imphcated in RDEB, before 
linkage analyses excluded its causative role (Colombi et aU 
1992). In the present study, we found overrepresentation of 
MMPl and MMP3 in CM of RDEB cells, and of MMP2 and 
MMP14 in the ECM. This was surprising, since increased 
proteolytic activity would apparently contrast with the 
accumulation of ECM in RDEB skin. The contradiction can 
be explained, at least in part, by the observation that although 
MMP 14 protein abundance was increased, protease activity 
was significantly reduced in RDEB cells. In vitro, the RDEB 
cells mimicked the phenotype of MMPi 4-deficient mouse 
embryonic fibroblasts, i.e., both cells were unable to degrade 
collagen I gels (Holmbeck et al, 1999). In RDEB, decreased 
proteolysis could be caused by increased inhibition or 
processing of MMP14, by hampered turnover, or a combina- 
tion of events. In support of the first prediction, increased 
levels of TIMP3 were identified in RDEB ECMs, and the 
addition of recombinant TIMP3 to control samples blocked 
collagen I degradation. In contrast, the addition of RDEB CM, 
which contained soluble protease inhibitors, to control cells 
had no effect rather pointing to a matrix-associated inhibitor 
like TIMP3 (Leco et al, 1994). The second prediction is 
supported by the fact that the ectodomain of MMP14 is shed in 
RDEB cells, rendering the membrane-bound stump inactive. In 
addition, MMP14 protein turnover was altered. Therefore, 
reduced MMP14 activity is a result of a combination of events. 

Apart from the fact that the three points above will open new 
avenues of research into disease mechanisms in RDEB, a very 
exciting observation in this first disease proteomics study on 
RDEB was that the proteome profiles correlated with the 
severity of RDEB during disease progress. Specifically, ECM 
proteome alterations of patient 2, who developed the most 
severe RDEB phenotype in the course of the disease, were most 
different from those of the other patients. The clinical 
manifestations were intermediate in case of patients 1 and 3 
and the respective ECM proteomes were most similar; the 
mildest clinical manifestations of the disease were observed in 
patient 4. These intriguing findings hold promise for future 
apphcations in personahsed medicine, e.g., by developing 
prognostic markers based on the classification of protein 
abundances. Culture of skin fibroblasts is an easy and 
minimally invasive procedure, and quantitative proteomics 
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is a powerful technique to identify cellular features, which are 
pertinent for both prognostication and development of 
targeted therapy approaches on an individual basis. 

Taken together, the present study emphasises the potential 
of disease proteomics approaches to identify novel mechan- 
isms of genetic diseases. The quantitative protein abundance 
information of extracellular proteins of healthy and patholo- 
gically altered cells provides a rich information source for the 
scientific community highlighting the complex changes in 
cellular microenvironment due to loss of a single protein, C7. 
The proteome alterations in RDEB serve as a starting point for 
discovery of hitherto unknown molecular mechanisms that 
not only explain the clinical phenotypic variation, but also 
identify targets for therapeutic intervention. 

Materials and methods 
Patients and RDEB diagnosis 

Skin specimens were obtained from four newborns with skin 
blistering, but a negative family history for skin diseases. Immuno- 
fluorescence mapping revealed lack of C7 in the skin in all patients. 
Mutation analysis of the C7 gene, C0L7A1 (Kern et al, 2006, 2009) 
disclosed mutations leading to premature stop codons and thus 
confirmed the diagnosis (Table I) . The clinical follow-up of the patients 
was during 2.5-10 years at the Epidermolysis Bullosa Center Freiburg. 
The study was approved by the Ethics Committee of Freiburg 
University and conducted according to the Declaration of Helsinki. 



Cell culture and stable isotope labelling in cell 
culture 

Primary dermal fibroblasts were isolated from the skin of the four 
patients and, as controls, from the foreskin of circumcised 3-, 4- and 9- 
year-old healthy boys (Sprenger et al, 2013). The cells were 
subcultured and passaged in SILAC-DMEM (Thermo Fisher, Langen- 
selbold, Germany), supplemented with 10% dialysed FBS (Gibco, 
Invitrogen), 1% Antibiotic-Antimycotic (100 x , Invitrogen), 1% l- 
glutamine (PAN Biotech), 42 mg/1 L-arginine (Sigma- Aldrich), 73 mg/ 
1 L-lysine (Sigma-Aldrich) and 82 mg/1 proline (Sigma- Aldrich) for 
the control population. RDEB fibroblasts were cultured and fully 
labelled for 2 weeks in the same way but through incorporation of 
either L-arginine-^^C6'*N4 and L-lysine-^H4 (Arg6, Lys4; Sigma-Aldrich) 
or L-arginine-^^C6-^^N4 and L-lysine-^^C6^N2 (ArglO, Lys8; Sigma- 
Aldrich). 



CM collection and ECM isolation 

In all, 5 X 10^ fully SILAC-labelled fibroblasts were seeded on 10 cm^ 
cell culture dishes (BD, Heidelberg, Germany). One day later, the 
medium was changed and 50|ig/ml ascorbate (Sigma-Aldrich) was 
added to allow proper folding and secretion of collagens (Myllyharju 
and Kivirikko, 2004). Before harvesting the CM, the cells were washed 
3 X with PBS and starved for 14 h in serum-free medium to avoid 
interference of serum proteins with the MS analysis. CM was collected 
on day 4, centrifuged at 500g 4°C for 10 min and filtered using 0.2 |im 
filters to remove cells and debris to prevent contamination by 
intracellular proteins. The CM of control and RDEB fibroblasts were 
mixed in a 1:1:1 ratio (the ratio was based on cell number) and 
concentrated by ultrafiltration using vivaspin columns with an mwco 
of 10 000 Da. Approximately 150 jig total protein (50 |ig of each sample) 
was loaded on an SDS gel. 

After 2 days in DMEM containing 50 |ig/ml ascorbate and 10% FCS, 
the fibroblasts were washed three times with PBS. The cells were 
removed from the underlying ECM with 0.5% Triton X-100 in 20 mM 
NH4OH for 30 s (Vlodavsky, 2001). The remaining ECM was 
then carefully washed 5 x with PBS to eliminate intracellular 



contaminants. The ECM was solubihsed with 4% SDS in 0.1 M 
Tris-HCl, pH 7.6 (Wisniewski et al, 2009). Aliquots of the different 
ECM isolations were mixed 1:1:1 (~50|ig for each sample) and 
prepared for MS analysis to determine the accurate mixing ratios. 
Based on this, whole ECM samples were mixed and concentrated using 
vivaspin columns with an mwco of 10 000 Da. 



MS sample preparation 

Samples were lysed in SDS-PAGE loading buffer, reduced with 1 mM 
DTT (Sigma-Aldrich) for 5 min at 95°C and alkylated using 5.5 mM 
iodoacetamide (Sigma-Aldrich) for 30 min at 25°C. Protein mixtures 
were separated by SDS-PAGE using 4-12 % Bis-Tris mini gradient gels 
(NuPAGE, Invitrogen). The gel lanes were cut into 10 equal slices, 
which were in-gel digested with trypsin (Promega, Mannheim, 
Germany) (Shevchenko et al, 2006) , and the resulting peptide mixtures 
were processed on STAGE tipps as described (Rappsilber et al, 2007). 
In case of the protomap approach, gel lanes were cut into 17 equal 
shces and processed as mentioned above (Dix et al, 2008). 



Mass Spectrometry 

Mass spectrometric measurements were performed on LTQ Orbitrap 
XL mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) 
coupled to an Agilent 1200 nanoflow-HPLC (Agilent Technologies 
GmbH, Waldbronn, Germany) . HPLC-column tips (fused silica) with 
75 |im inner diameter (New Objective, Woburn, MA, USA) were self 
packed (Gruhler ef a/, 2005) with Reprosil-Pur 120 ODS-3 (Dr. Maisch, 
Ammerbuch, Germany) to a length of 20 cm. Samples were applied 
directly onto the column without pre-column. A gradient of A (0.5% 
acetic acid (high purity, LGC Promochem, Wesel, Germany) in water 
(HPLC gradient grade, Mallinckrodt Baker B.V., Deventer, The 
Netherlands) and B (0.5% acetic acid in 80% ACN (LC-MS grade, 
Wako, Germany) in water) with increasing organic proportion was 
used for peptide separation (loading of sample with 2% B; separation 
ramp: from 10-30% B within 80 min). The flow rate was 250nl/min 
and for sample application 500nl/min. The mass spectrometer was 
operated in the data-dependent mode and switched automatically 
between MS (max. of 1 xlO^ ions) and MS/MS. Each MS scan was 
followed by a maximum of five MS/MS scans in the linear ion trap 
using normahsed collision energy of 35% and a target value of 5000. 
Parent ions with a charge state fromz = 1 and unassigned charge states 
were excluded for fragmentation. The mass range for MS was 
m/z = 370-2000. The resolution was set to 60 000. Mass-spectrometric 
parameters were as follows: spray voltage 2.3 kV; no sheath and 
auxiliary gas flow; ion-transfer tube temperature 125°C. 

Identification of proteins and protein ratio 
assignment using MaxQuant 

MS raw data files and respective information can be downloaded from 
the Peptide Atlas Data Repository (http://www.peptideatlas.org/ 
PASS/PASS00080) (Desiere et al, 2006). 

The MS raw data files were uploaded into the MaxQuant software 
version 1.2.0.18 (Cox and Mann, 2008), which performs peak 
detection, SILAC-pair detection, generates peak lists of mass error 
corrected peptides and database searches. A full-length IPI human 
database containing common contaminants such as keratins and 
enzymes used for in-gel digestion (based on IPI human version 3.68) 
was employed, carbamidomethylcysteine was set as fixed modifica- 
tion, oxidation of methionine, proline and lysine, and protein amino- 
terminal acetylation were set as variable modifications. Triple SILAC 
was chosen as quantitation mode. Three miss cleavages were allowed, 
enzyme specificity was trypsin/P + DP, and the MS/MS tolerance was 
set to 0.5 Da. The average mass precision of identified peptides was in 
general < 1 p.p.m. after recalibration. Peptide lists were further used 
by MaxQuant to identify and relatively quantify proteins using the 
following parameters: peptide and protein FDRs, based on a forward- 
reverse database, were set to 0.01, maximum peptide posterior error 
probability (PEP) was set to 0.1, minimum peptide length was set to 6, 
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minimum number peptides for identification and quantitation of 
proteins was set to two of which one must be unique, minimum ratio 
count was set to two, and identified proteins have been requantified. 
The 'match-between-run' option (2min} was used. In case of the 
protomap approach, spectral counts were used as an estimate for 
relative protein quantification. 



Data analysis 

To obtain a list of potential secreted/extracellular proteins, the data 
were filtered using the freely available Perseus software (Cox et al, 
2011) to extract proteins with the already annotated GO terms: 
extracellular region, extracellular space and ECM. Proteins that were 
present in the ECM fraction and CM were also included in the final 
protein lists, as well as proteins predicted by in silico analyses to carry 
a signal peptide (SignalP; default settings, www.cbs.dtu.dk/services/ 
SignalP/) (Petersen et al, 2011). Subsequently, data were readjusted to 
a median of 1 and log2 transformed. DAVID 6.7 BETA (Huang da et al, 
2009) was used to investigate protein enrichment in the ECM and CM 
fraction using functional annotation clustering with respective SP and 
PIR keywords with the default settings with a minimum significance of 
P<0.01. 

The seven fibroblasts samples (pts and ctrls) were hierarchically 
clustered by the abundance of the identified proteins and the ratios for 
the individual proteins were k-means clustered to partition the 
proteins to six discrete clusters. To address the biological implications 
of the proteins in each cluster. Biological Process and Molecular 
Function GO terms were retrieved using DAVID (Huang da et al, 2009) 
and binomial probability test followed by Benjamini and Hochberg 
P-value adjustment was used to extract terms enriched in each cluster 
when tested against the remaining five clusters. A P-value after 
adjustment below 0.05 and at least five occurrences in the cluster were 
required to regard the enrichment as significant. 

Received hydroxyproline and hydroxylysine sites were extracted for 
high confidence sites by requiring a MaxQuant localisation probability 
of >0.75 and an Andromeda Score of ^100 (Olsen et al, 2006; Cox 
et al, 2011) . Pre-aligned peptides were further analysed by motif-x vl .2 
10.05.06 (Schwartz and Gygi, 2005; Chou and Schwartz, 2011) with a 
central character of either K for hydroxylysine or P for hydroxyproline 
sites, a width of 13, a series of occurrences of 10 and a significance of 
0.000001. 



Western blot 

Confluent cell cultures were extracted with 1% Nonidet P-40, 0.1 M 
NaCl, 0.025 M Tris/HCl, pH 7.4, containing Complete Protease 
Inhibitor Cocktail (Roche). The extracts were cleared from nuclei by 
centrifugation at 16 000 g for 10 min at 4°C, and the protein content of 
the extracts was determined by DC Protein Assay (Bio-Rad, Miinchen, 
Germany). An ahquot containing 50|ig of protein was separated on 
SDS-PAGE using a 7.5-12.5% gradient gel and transferred onto a 
nitrocellulose membrane. The membranes were blocked with 5 % milk 
powder in 1 x TBS with 0.1 % Tween-20 for a minimum of 1 h at room 
temperature and incubated for 1 h or overnight with one of the 
following antibodies: anti-nidogen (sc-33141), anti-fibulin-1 (sc- 
25281), anti-laminin p-1 (sc-5583), anti-laminin y-1 (sc-17751), anti- 
actin (sc-47778) all from Santa Cruz Biotechnology; anti-Tenascin-C 
(MAB2138; R&D Systems, Wiesbaden, Germany), anti-MMP-14 
(ab51074; Abeam), anti-Big-h3 (D31B8; Cell Signaling Technology), 
or the collagen VII antibody NC2-10 (Bruckner-Tuderman et al, 1995). 
HRP-conjugated secondary antibodies and a chemiluminiscent detec- 
tion assay (Immobilon Western, Millipore, Schwalbach, Germany) 
were used for visualisation according to manufacturer's instructions. 



Immunofluorescence staining 

For indirect immunofluorescence staining, fibroblasts grown on 
coverslips or skin cryosections were fixed with 4% PFA, or cold 
acetone and blocked with 2% BSA in PBS for 30 min at room 
temperature, followed by incubation with the primary antibody diluted 
in 0.2% BSA in PBS overnight at 4°C. After incubation with the 



secondary antibody in PBS for 1 h, the samples were embedded in 
fluorescence mounting medium (Dako, Hamburg, Germany). For 
staining of the ECM, 5 x 10"^ fibroblasts were grown on coverslips in 
DMEM containing 50 |ig/ml ascorbate for 1 week. The cells were then 
removed as described above, and the ECM was fixed with 4% 
paraformaldehyde in PBS for 15 min and permeabihsed with 0.1% 
Triton in PBS for 5 min at room temperature. After blocking for 30 min 
in 2 % BSA in PBS, the ECM was incubated with the primary antibodies 
diluted in 1 % BSA in PBS for 1 h. After incubation with the secondary 
antibody, the ECM was washed in PBS and embedded in fluorescence 
mounting medium (Dako). The following primary antibodies were 
used: anti-fibronectin (AB26245, AB2413; Abeam), anti-tenascin-C 
(MAB2138; R&D Systems), anti-Big-h3 (D31B8; Cell Signaling Technol- 
ogy), anti-nidogen (sc-33141), anti-Laminin-p/yl (Progen, 10765), 
anti-ColIV (AB15633, Abeam) and the collagen VII LH7,2 (Bruckner- 
Tuderman et al, 1995). Alexa488- or Alexa594-conjugated secondary 
antibodies (Invitrogen) were used. Pictures were taken with an IF 
microscope (Zeiss Axio Imager; Zeiss, Oberkochen, Germany) . 



Staining of paraffin-embedded skin biopsy 
specimens 

Skin specimens from the extremities of five RDEB patients (collagen 
VII negative, aged 25-45 at time of sampling) and five healthy controls 
(age and localisation matched) were Elastica van Gieson stained using 
standard methods. Mean dermal staining intensity in grey-scale 
pictures (16 bit) was measured with Image J (http://rsbweb.nih. 
gov/ij/) and normalised to mean epidermal staining. Two-tailed, 
unpaired T-test was used to calculate significance. 

In situ collagen degradation assay and 
zymography 

In all, 12- well culture plates were coated with a layer of neutralised rat 
tail collagen I (Gibco, Invitrogen) with a concentration of 300|ig/ml 
(1 ml/well) as described by the supplier and incubated for 2 h at 37°C 
to allow fibril formation. After washing with water and drying of the 
collagen gels, 5 x 10"^ fibroblasts were seeded onto the film. Five hours 
later, culture medium (DMEM, 10% FCS, 1% l-G1u, 1% penicillin/ 
streptomycin) was added. The following day, culture medium was 
changed to serum-free DMEM and TNF-a, ConA or protease inhibitors 
were added. Following concentrations were used: lOng/ml TNF-a 
(Peprotech, Hamburg, Germany), 2nM ConA (Peprotech), 25|iM 
TAPI-0 (CalBiochem), 1 mM PMSF (Sigma- Aldrich) , 1 |iM SBTI 
(Sigma- Aldrich) , 1 |ig/ml Leupeptin (Sigma- Aldrich) , 50|iM Furin 
Inhibitor I (Decanoyl-Arg-Val-Lys-Arg-CMK, Calbiochem/Merck, 
Darmstadt, Germany), 100 nM human recombinant TIMP3 (R&D 
Systems). After an incubation of 3 days in 37°C CM was harvested 
and the cells were removed by using 0.05% Trypsin and 0.02% EDTA 
in PBS (PAN Biotech, Aidenbach, Germany) . The remaining collagen 
was visualised by Coomassie Brilliant Blue R-250 (Bio-Rad). Images 
were captured with an Image Quant™ LAS 4000 mini (GE Healthcare, 
Freiburg, Germany). 

For MMP2 zymography, CM was centrifuged at 500 g at 4°C for 
10 min and cells and cell debris of the remaining supernatant were 
removed by filtration using 0.2 |im filters. In all, 5 |ig of proteins was 
separated on 10 % ready gel gelatin zymogram gels (Bio-Rad) for 3 h at 
100 V without prior reduction. Subsequently, the gel was washed 4 x 
in 2.5% Triton X-100 to remove SDS and incubated overnight in 
zymogram developing buffer (Bio-Rad) at 37°C to allow gelatin 
digestion. Gelatinase activity was visualised as colorless bands after 
staining the gel with Coomassie Brilliant Blue R-250. 



siRNA-mediated silencing of collagen VII 
expression 

For collagen VII knockdown, primary human dermal fibroblast were 
seeded at 30-50% confluence and transfected with siRNA duplexes 
specific to collagen VII the following day (Martins et al, 2009) or with 
irrelevant control siRNA duplexes (Eurogentec, Liege, Belgium) and 
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Lipofectamine 2000 (Invitrogen) according to manufacturer's recom- 
mendations. Collagen VII expression was analysed by western blot. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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